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ABSTRACT 

Online Digital Sky Survey (DSS) material is often used to obtain information on newly 
discovered variable stars for older epochs (e.g. Nova progenitors, flare stars, . . . ). We 
present here the results of an investigation of photometry on online digital sky survey 
material in small fields calibrated by CCD sequences. We compared different source 
extraction mechanisms and found, that even down near to the sensitivity limit, despite 
the H-compression used for the online material, photometry with an accuracy better 
than 0^1 rms is possible on DSS-II. Our investigation shows that the accuracy de- 
pends strongly on the source extraction method. The SuperCOSMOS scans, although 
retrieved with an higher spatial resolution, do not give us better results. The methods 
and parameters presented here, allow the user to obtain good plate photometry in 
small fields down to the Schmidt plate survey limits with a few bright CCD calibra- 
tors, which may be calibrated with amateur size telescopes. Especially for the events 
mentioned above, new field photometry for calibration purposes mostly exists, but 
the progenitors were not measured photometrically before. Also the follow up whether 
stellar concentrations are newly detected clusters or similar work may be done with- 
out using mid size telescopes. The calibration presented here is a "local" one for small 
fields. We show that this method presented here gives higher accuracies than "global" 
calibrations of surveys (e.g. GSC-II, SuperCOSMOS and USNO-B). 

Key words: methods: data analysis - techniques: photometric - surveys. 



1 INTRODUCTION 

The photographic sky surveys produced by large Schmidt 
telescopes have proved to be among the most useful and 
enduring weapons in the astronomical armory. With the 
start of the Digitized Sky Survey (hereafter DSS; Lasker 
& McLean 1994) new powerful tools have been estab- 
lished. To provide convenient access to these data, the im- 
ages have been compressed using a technique based on the 
H-transform to reduce the data volume. Although the tech- 
nique is lossy, it is adaptive so that it preserves the signal 
very well. They typically compress the data by a factor of 7, 
but much higher compression ratios are possible. The south- 
ern surveys are covered by the SuperCOSMOS Sky Survey 
(hereafter SSS; Hambly et al. 2001a), too. This survey pro- 
vides us with images at even higher spatial resolution. While 
there are extensive studies on the photometric calibration of 
whole plates or even surveys on original uncompressed data 
(e.g. Reid & Gilmore 1982, Lasker et al. 1990, Russell et al. 
1990, Morrison et al. 1997, Monet et al. 2003 = USNO-B), 
the possibilities of local small field calibration on compressed 
online available material were not detailed studied for stel- 
lar photometry. Galaxy surveys use this material more often 
and thus much better studies are available there (e.g. Kron 
1980, Spagna et al. 1996, Hambly et al. 2001b). Although 



there is often the need to look up the older epoch data, 
especially for eruptive variable stars, only limited effort on 
the source extraction and calibration technique for stars was 
done (e.g. King et al. 1981, Humphreys et al. 1991, Hortnagl 
et al. 1992, Kimeswenger & Weinberger 2001, Andersen & 
Kimeswenger 2001, Kimeswenger et al. 2002a, 2002b, 2003). 
The goal of this study here is to show how to optimize 
source extraction parameters in small fields around inter- 
esting sources. This allows a user to derive a high quality 
local calibration of online digitized plates, having only a few 
bright calibration sources in the field. Such calibrators can 
be obtained easily by CCD cameras at amateur size tele- 
scopes or training equipment of universities (Kimeswenger 
2001, Bacher et al. 2001, Lederle & Kimeswenger 2003). 



2 DATA 

The photometric data was obtained from the online servers 
at the STScI ( http:/ /archive. stsci.edu/cgi-bin/dss.plate^finder^ , 
ESO {http://archive.eso.org/dss/dss) and ROE 

{http://www-wfau.roe.ac.uk/sss/pixel.htnil). The charac- 
teristics of these surveys are summarized in Table Q 

As photometric sequences those of Henden (IMT^ 
were chosen. These widely used sequences were compared 
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Table 1. Parameters and colour equations for sky survey plates used here. Where several colour equations are given, always the first 



was used here ( 


see text). Some SuperCOSMOS scans are 


only partly available yet 




survey 


plate/filter 


band 


DSS scan 


SuperCOSMOS 


colour equation 


reference 








resolution 


resolution 






Pal-QV No 


IIaD+W12 


V 


V/7 








SERC-J/EJ 


IIIaJ+GG395 


Bj 


V/7 


0'/67 


Bj = B - 0.28 (B-V) 


Blair & Gilmore 












Bj = B - 0.20 (B-V) 


King et al. IT?^ 












Bj = B - 0.23 (B-V) 


Kron frM)l 


POSS II J 


IIIaJ+GG385 


Bj 


V/0 




Bj = B - 0.28 (B-V) 


Reid et al. IT?M1 


POSS I E 


103aE 


E 


r/7 


0'/67 


E = Rc 


Spagna et al. |1996t 


SERC-ER/R 


IIIaF+OG590 


r = R59 


V/0 


0'/67 


r = Rc + 0.04 (V - Rc) 


Hortnagl et al. IllMI 


POSS II R 


IIIaF+RG610 


r = Rei 


V/0 




used same as SERC-R 





and tested with own measurements from the Innsbruck 
60cm telescope (Kimeswenger et al. 2002a, Lederle & 
Kimeswenger 2003) around V838 Mon and CI Aql and with 
those by K.S. from ESQ NTT around V4332 Sgr. 

Colour equations for plate/filter combinations were ob- 
tained by various authors in different manner. While King 
et al. (|1981|) have arbitrarily chosen a relation based on vi- 
sual estimate from the passbands, Hortnagl et al. ( 1992 ) and 
Henden IFM)^ for use in USNO B (Monet et al. 2003) used 
spectrophotometric catalogues to fold with published filter 
+ detector response curves. Kron STM)\. Blair & Gilmore 
(UnHS, Reid et al. (fTMTll . Humphreys et al. ITMT)l and 
Spagna et al. ("1996") used photoelectrical or CCD photome- 
try of stars in standard bands to derive the equations. Where 
available we used the latter ones here. We kept the native 
magnitude system of the plates and converted the CCD mag- 
nitudes to this system since the POSS-II and the southern 
plates for a given field were taken at different epochs, so vari- 
ables could get really weird colours, along with some other 
factors. 



\ ' \ ' \ ' r 
" * DSS II red (plate ER772)- fixed aperture 5" 
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Figure 1. A sample of photographic (DSS-II red) fixed aperture 
(5^^) magnitudes vs. CCD magnitudes. While the faint end gives 
us a linear relation, the core saturation of the bright stars lead to 
an asymptotic convergence to a constant value at the bright end. 



3 CALIBRATION METHOD 

The source extraction was performed using SExtractor V2 
(Bertin & Arnouts 1996). We used aperture photometry 
with fixed sized apertures and Kron first moment adaptive 
radius apertures (Kron 1980) as implemented to SExtractor. 
Due to the core saturation of bright stars we expect a curve 
reaching asymptotically a constant value at the bright end 
and being quasi linear at the faint end (see Figure 0. 

Such a curve was " invented" for the calibration of pho- 
tographic densities by Moffat (1969 ). He investigated in de- 
tail the calibration curves and the image growth of the stars 
on Hamburg Schmidt plates. We here use the equation of 
the same form 

mphot = ao X log (iO"ix('"ccD-a2) + i) + ^3 (1) 

where mphot is the magnitude coming directly from the 
source extraction on the DSS or SuperCOSMOS plate, 
mccD is the CCD magnitude corrected to the colour system 
of the plate and ai; i ^ [0, 3] are free fit parameters. As we 
measure in case of stars an integrated value 

/"-Raperture /"-Raperture 

/ /(T(r))rdr- / /(S(r))rdr (2) 

Jo Jo 

were T(r) is the really transmitted light in the measurement 



machine for the star and S(r) that of the sky background. In 
the software we use the sky background is calculated as 2D 
polynomial in the region around the stellar aperture after 
exclusion of contamination by other sources. Details of this 
procedure can be found in Bertin & Arnaults (1996). This 
transmission is a function of seeing - and as shown by Moffat 
(1969) and Kimeswenger (1990) - strongly changes due to 
the wavelength dependent scattering within the plate. As 
shown already by Moffat and more in detail by Irwin & Hall 
(1983) simple fits of this growth function with e.g. Gaussians 
(as working fine for direct imaging CCDs) do not work well. 

Assuming ai ^ 1 leads to the linearized "Baker den- 
sities". This works fine in case of some emulsion and wave- 
length domains (Hoffman et al. 1998). This fact was widely 
used to "linearize" the calibration curves (e.g. the COSMAG 
values in the SuperCOSMOS catalogue) and to correct for 
the remaining deviation by semi-empirical functions (Ham- 
bly et al. 2001a). In fact this leads to "global" parameters 
for the liberalization and free parameters fitting deviations. 
This is very suitable for global calibrations but not needed 
for the intended local calibration here. 

But our main uncertainty in case of the public DSS 
data in fact is the unknown conversion function / in the 
integration (Equ. 2) above. We measured the value of the 
unexposed plate (thereafter called chemical fog) and found 
strong variations from plate to plate. The machines were 
not calibrated in the same way all the time when scanning 
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DSS2. Thus a general clear correlation of our fit parameters 
with the theoretical values by Moffat (1969) is not possible. 
This seems to be a major difference to the SuperCOSMOS 
scans, were all plates seemed to be scanned and calibrated 
in the same way. Therefore there is no way around a fitting 
of the parameters to obtain a proper solution. 

While as gives an estimate for the core saturation value, 
which strongly depends on the aperture size, ao x ai rep- 
resents the inclination of the linear section, ai also influ- 
ences the curvature around the knee and a2 gives more 
or less the start of the linear section. An individual fit 
^phot(^ccD; ao, ai, a2, as) for each CCD sequence and each 
aperture was derived. The correlation coefficient for all fits 
was always > 0.995. We derived the standard deviation 
of the fitted values af in bins of Am = 0™1 for m > 12"'0. It 
is nearly negligible for bright sources and increases steadily 
towards the exposure limit. Near the limit the error will orig- 
inate partly also from the CCD sequences (see Figure |2j. 

The reverse function of (1) then has to be used to derive 
the mccD of targets 



mccD = — X log 
ai 



10" 



(^phot ~°-3) 



1 +a2 



(3) 



To derive an error estimate for the resulting magnitudes an 
for the reverse function we used the derivative 



/ dmphot \ " 
VdmccD / 



aoai 



(l + lO""^^ 



(4) 



This clearly leads to a very high effective error at the 
bright end of the calibration where the saturation effects the 
data. There the isophotal magnitudes are definitely better. 
But as we have shown in the introduction we are not in- 
terested in this part of the calibration as it can be reached 
better by other means. In Figure |21the rms residuals for a fit 
are shown. For most of the brightness region of interest the 
error is about 0™1 or even below. Finally we tested whether 
the fits, respectively the residua are a function of the stellar 
colour. This test gives us information on the quality of the 
colour equations used. 

The smaller the aperture we used, the smaller the error 
was (Figure El at higher magnitudes. On the other hand 
the saturation goes deeper. Often adaptive Kron apertures 
(Kron 1980) are used to avoid saturation. This method is 
also implemented in SExtractor. Comparing our result to 
the accuracy shows that the rms is larger at all wavelengths. 
Crowding seems to confuse too often the adaptive apertures. 
Kron adaptive apertures, using two or more threshold levels, 
are indeed a kind of profile fitting. As already pointed out 
by Irwin & Hall (1983), profile fitting leads to the highest 
errors. They find best results for isophotal photometry, while 
we find in our crowded fields in the galactic plane best results 
by using small apertures. 

The solutions for the individual scans strongly differ 
from plate to plate. There is no general solution for the pa- 
rameters. As shown in Figure 0] the plates are exposed and 
scanned so differently that neither the saturation (as) not 
the inclination at the higher magnitudes (ao x ai) nor the 
curvature (ai) or the position of the knee (a2) is fixed some- 
how. Normally one does not have deep CCD sequences as 
we used here for the base study. Thus we tested the ro- 
bustness of the fits with a small number of bright CCD 
calibrators. But in fact we were searching for an extrapo- 



0.2 



0.1 




13 14 15 16 17 18 19 
m(R) 

Figure 2. The rms as function of magnitude for plate ER772 in 
the region around the variable star V838 Mon. The error mostly 
lies at or below 0™1. The slow increase towards the fainter end is 
partly denoted to the accuracy of the photoelectric sequence of 
Henden (2002). At the bright end the error increases rapidly due 
to the core saturation. The smaller "jump" (at 14™5) is due to a 
"cluster" of slightly blended stars around the target V838 Mon, 
having all about the same apparent magnitude. They also might 
have inaccurate CCD magnitudes, as the calibration sequence was 
taken during a bright outburst of the Target. 



lation if only a few bright CCD calibrators are available. We 
thus tried to derive some of the parameters by basic mea- 
surement on the plate itself. For this purpose we obtained 
the value of the unexposed plate If at the plate edge (often 
called chemical fog) , the value of the sky background Is and 
finally that of the overexposed plate lo in the center of a 
bright star nearby. All those measurements are used " as is" 
in the published files. These units are the original natural 
analog/digital converter units folded by the (to the user un- 
known) function / (see equation 2). Looking to equation (1) 
one can see that as gives the magnitude of a star were the 
overexposed core is just as big as the (fixed size) aperture. 
Inserting the saturation and the sky level to equation (1) 
this leads together with the aperture d (in units of pixels) 
to: 



as 



-2.5 log 



{lo - Is) 



(5) 



In our tests we found a strict correlation of the chemical 
fog and the sky level with the parameter a2. The dynamic 
range from the chemical fog to the sky level {Is — If) gives 
us 



mphot(a2) - -2.5 log 



— {Is 



(6) 



The CCD magnitude of a star with mphot(a2) (or the mean 
of a set of stars within mphot(a2) ±0™05) is equal to a2. This 
gives a limit for the deepness of the CCD sequence needed 
for the calibration. 

Finally inserting the CCD magnitude of the curvature point 
a2 into equation (1) we get with the two equations (5) and 
(6) above a solution for ao 



ao = 



2.5 

"log(2) 



log 



Is -If 



lo-Isi 



(7) 



This parameter ao is independent from the aperture, as it 
describes the inclination of the curve at stars much smaller 
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Figure 3. The deviation of fits with different apertures in the 
15'xl5'field around V838 Mon. 
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Figure 4. The field around the variable star AM Her on the plates 
XJ229 and XJ230. The overlapping region allows us to study the 
different behaviour of the plates individually. 



than the aperture. As we now were able to fix 3 (out of 4) 
parameters, the fit of the remaining parameter ai was robust 
against reducing the number and the limiting magnitude of 
the calibrators. 

Figure |5] shows that the limiting magnitude for the CCD 
sequence has to be deep enough to reach the value of a2 
derived above. Deeper sequences do not improve the quality 
of the parameter anymore. Thus a real extrapolation from 
bright CCD sequences to faint stars is possible (see Figure 
6). 



0.5- 



^ 0.4- 



0.3- 



4 THE H-COMPRESSION 

The DSS data is compressed by the use of a "lossy" algo- 
rithm. The H-compression (White et al. 1992) uses typical 
factors of 5 in crowded fields and up to 10 at high galac- 
tic latitudes (Lasker et al. 1996). We investigated the in- 
fluence of the H-compression in one field around the sus- 
pected symbiotic variable star V471 Per. For comparison 
the plate from the SuperCOSMOS data archive, scanned at 
much higher resolution of 0^67, was median filtered 3x3 pix- 
els (= lV98xlV98) and then rebinned to the DSS pixel size 



I I I I I I I I I I I 

10 12 14 16 18 20 

^CCD 

Figure 5. The evolution of the fits of parameter ai as func- 
tion of the limiting magnitude. It converges to the optimum at 
magnitudes around a2 which was derived independently before as 
described above. 
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Figure 6. The real 4 parameter fit is nearly identical with the fit 
fixing 3 parameters and using the sequence up to mccD = 15™ 0. 
Even the fit with very bright calibrators only, although systemat- 
ically wrong at the faint end, does not deviate more than about 
0™5 - the rms found often for globally calibrated astrometric sur- 
veys like e.g. USNO B. 




14 16 18 20 

^from H-compressed DSS2 



16 18 20 

from H-compressed DSS2 



Figure 7. The photometry achieved on uncompressed (Su- 
perCOSMOS plate rebinned to DSS resolution) vs. DSS-II H- 
compressed images around the suspected variable star V471 Per 
using the method described in section 3. There is obviously nei- 
ther a systematic effect nor any degrading of the photometric 
quality. 



of 1^/01. A CCD calibration sequence is given by Henden & 
Munari (|2UU1|1 . The result in Figure Q shows that there is 
no degeneration of the photometry even down to the plate 
limit. The rms of the regression is below 0™09 at m < 20™ 
and increases to 0™14 at 20 ^ m < 21™5. 



5 COMPARING TO SURVEY CALIBRATIONS 

Recently a set of globally calibrated sur- 
veys was released. While the GSC-II 
\http://www-gss's. stsci.edu/gsc/gsc2/ GSC2home.htm ) is 
based on the same scans (but before H-compression) we 
used for our investigation, the USNO-B catalog (Monet et 
al. 2003) and the SSS catalog (Hambly et al. 2001a) used 
data from different scanning machines. The USNO-B only 
uses an 8-bit linear converter for the transmitted light. 
This undersamples the greyscale variations in the centre 
of the stellar images. As already pointed out by Monet et 
al. (2003), the photometry never was a main goal of this 
survey. The SSS uses scans with a linear 15-bit greyscale 
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Figure 8. The deviation of the magnitudes given in the USNO-B 
catalog (Monet et al. 2003), the GSC 2.2 (STScI, 2001), and the 
SSS catalog (Hambly et al. 2001a) from those of Henden (2002) 
in the 15^xl5^field around AM Her. 



from a CCD. The data is pre-calibrated using an average 
slope for the so called linear part of the emulsion of 7 = 2. 
Then a complex post-processing including CCD calibrators 
and colors is started (Hambly et al. 2001a). This is essential 
for wide area statistical studies, but not needed in our local 
calibration of individual sources. The comparison of the 
Henden (2002) CCD sequences with these catalogues for 
two of our fields is shown in Figs.|H]&|3| 

The GSC2.2 do show a much higher noise but (nearly) 
no systematic effects. The calibration is not well documented 
in the literature neither at the survey homepage. Manual in- 
spection of the sources in our fields leads us to the conclusion 
that at least a part of this error is originating from improper 
deblending of faint sources. The typically higher crowding 
on the red survey plates thus causes a higher error there. 

The USNO-B photometry is based on extrapolation of 
bright TYCHO stars and a few faint calibrators on some 
plates. Finally the remaining plates were adjusted using the 
overlaps. We find generally very high (^ 0™5) systematic 
deviations for all our fields. 

The SSS uses a grid of faint CCD calibrators fields to 
define an average gradient of the isophotal magnitudes vs. 
CCD magnitudes. These were used to extrapolate towards 
the plate limit the individual field- by- field calibration curves 
derived from bright standards. Post-processing by a com- 
plex color classification and the plate overlaps were used to 
do fine-tuning (Hambly et al 2001a). The comparison with 
our results in the small fields leads to a much lower scatter 
than the surveys mentioned before. But systematic effects 
are sometimes evident. This is already pointed out by Ham- 
bly et al.: "The isophotal scale is nonlinear at the level of > 
0.5 mag; furthermore, such nonlinearities are not repeatable 
from field to field. " 
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Figure 9. The deviation of the magnitudes given in the USNO-B 
catalog (Monet et al. 2003), the GSC 2.2 (STScI, 2001), and the 
SSS catalog (Hambly et al. 2001a) from those of Henden (2002) 
in the 15^xl5^field around V838 Mon. 
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Figure 10. Left: The SSS avoids in the 'full-quality' catalogue 
the region around bright stars (shaded area). This causes targets 
like here V4332 Sgr being missed; Right: SSS and GSC-II do not 
deblend targets, which are properly done by the neuronal network 
software SExtractor (elliptical areas). Also faint stars having a 
good S/N ratio nearby bright stars are ignored (polygonal area). 



In a recent investigation of a field around the eruptive 
variable V4332 Sgr (Kimeswenger, 2005) another effect was 
crucial. The SSS, using the selection flags for high quality 
(= default), misses targets near bright stars (Figure QSJ- 
Also deblending, which works fine in SExtractor V2 (Bertin 
& Arnouts 1996), is not implemented to GSC-II and SSS 
properly. 

The calibration of the stars around V4332 Sgr by using 
ESQ NTT ODD images (limiting mag ^ 24?'0) obtained 
by SK gives, after manual cleaning blends and all stars in 
the vicinity of 10" around bright stars a very good result 
for GSC-II. Again the SSS has strong systematic deviations 
(Figure [TTJ. 
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Figure 11. The GSC-II and the SSS red magnitudes in a field 
around the eruptive variable V4332 Sgr (Kimeswenger, 2005). 
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Figure 12. Calibration of the field around V838 Mon using the 
CCD sequence (upper panel) and using GSC-II (lower panel). 



ALTERNATIVE CALIBRATION 
GSC-II 



CCD VS. 



As shown in the previous section the GSC-II, although hav- 
ing a higher rms error, shows nearly no systematic effects. 
Thus we tested a calibration by use of the GSC-II instead of 
the CCD sequence. This might be of interest in regions were 
no CCD calibrations are available or for a first quick result 
- e.g. especially after erruptive nova and nova like events. In 
Figure IT^ the results of a calibration using a CCD sequence 
and alternatively using only the GSC-II to derive the ai pa- 
rameters are shown. The results are promising. The system- 
atic effects are significantly below the noise. The accuracy 
of such a local "post-processing" of the GSC-II improves the 
photometry significantly. 



7 CONCLUSIONS 

The method we present here is well suited to calibrate lo- 
cally even H-compressed digital sky survey data by just us- 
ing a few bright calibrators. These calibrators easily can be 
achieved by using small amateur telescopes. The accuracy 
is about 0™1 for a wide range. Thus also colour-colour di- 
agrams are possible. This allows to do photometry of pro- 
genitors for eruptive events like novae as well as HR dia- 
grams when searching for new candidates for stellar clusters. 
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The importance for variable and eruptive star studies was 
pointed out recently by Munari et al. (j2002|l . Kimeswenger 
et al. (2002a) and Kimeswenger & Lechner ("2003"). Goran- 
skij et al. (2004) shows the strong implications on physi- 
cal models of such kind of DSS usage. As pointed out by 
Kimeswenger ("20051 several different epochs by using e.g. 
overlaps, not included all to the surveys, give us important 
information for such variable objects. Here we do not fo- 
cus on light curves of variables on a set of (similar) plates 
with respect to a reference plate. In such a case photomet- 
ric differential work is usually the best choice. Comparing 
single epoch data or data on different bands (e.g. POSS- 
I and POSS-II blue are different) with new CCD data the 
progenitors of eruptive variables indeed need absolute cali- 
bration. But also the investigations of small dark clouds by 
Stiiwe (1990 ), who used statistical calibrations with a BS 
model (Bahcall & Soneira 1984) and complementary galac- 
tic structure work like that of Curry & McKee (2000) and 
of Castellani et al. (|2()()l|l might be improved by this kind 
of caUbrations presented here. 

The local calibration presented here is clearly not overcom- 
ing global surveys (e.g. Hambly et al. 2001a or Monet et al. 
2003) for large statistical samples, but due to the higher ac- 
curacy (restricted to small fields) a perfect access for studies 
of individual targets is provided. Also the follow up whether 
stellar concentrations are newly detected clusters or similar 
work may be done without using mid size telescopes (e.g. 
Boeche et al. 2003). 
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